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ABSTRACT 


Large  single  crystals  of  benzene  (C5H5)  were  aligned  by  an 
optical  method.  The  aligned  crystals  were  pulsed  with  ultrasonic 
waves  at  a  frequency  of  approximately  10  megacycles  per  second  using 
quartz  crystal  transducers.  The  velocity  of  the  ultrasonic  waves 
was  measured  by  comparing  the  echo  transit  time  with  that  of  simul¬ 
taneously  generated  waves  moving  through  a  variable  length  path  in 
a  medium  in  which  the  velocity  of  sound  was  known  precisely. 

Velocities  of  longitudinal  and  transverse  waves  along  each  of 
the  principal  axes  of  the  orthorhombic  crystal  were  used  to  compute 
the  main  diagonal  terms  of  the  matrix.  The  off-diagonal  terms 
were  then  determined  by  velocity  measurements  along  arbitrarily 
selected  directions. 

Velocity  measurements  were  made  and  the  elastic  constants 
computed  over  the  temperature  range  170°K  to  250°K. 

Absorption  in  the  crystal  was  measured  over  the  same  temperature 
range  for  each  of  the  directions  used  for  velocity  measurements.  The 
absorption  was  determined  to  be  aelotropic. 

The  writers  wish  to  express  their  appreciation  for  the  assistance 
and  encouragement  given  them  by  Professor  0.  B.  Wilson  of  the  U.  S. 
Naval  Postgraduate  School  in  this  work. 
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INTRODUCTION 


Benzene,  the  fundamental  aromatic  hydrocarbon,  possesses  an 
orthorhombic  crystalline  structure  in  the  solid  state. 

The  low  melting  point  of  benzene  (5.2°C)  renders  the  use  of 
standard  diffraction  methods  for  the  determination  of  orientation 
difficult.  Benzene  crystals  are,  however,  transparent  to  visible 
light  and  possess  a  high  positive  birefringence;  thus  their  optical 
properties  provide  a  suitable  means  of  alignment  using  simple  optical 
methods . 

The  formulation  of  the  strain  energy  equation  for  an  orthorhombic 
crystal  results  in  nine  independent  elastic  constants(l) .  By  measure¬ 
ment  of  the  velocities  of  longitudinal  and  transverse  ultrasonic  waves 
along  the  three  principal  axes,  the  six  main  diagonal  terms  of  the 
matrix  can  be  computed.  The  remaining  three  elastic  constants  must 
then  be  determined  by  velocity  measurements  along  arbitrarily  selected 
directions,  using  the  previously  determined  values  for  the  main 
diagonal  terms . 

Using  the  techniques  for  alignment,  pulsing  and  velocity  measure¬ 
ment  described  in  this  thesis,  the  velocity  data  was  obtained  for  the 
computation  of  the  elastic  constants  over  the  temperature  range  170°K 
to  250°K. 

In  obtaining  the  velocity  data,  the  opportunity  was  taken  to 
observe  and  record  the  absorption  of  the  ultrasonic  pulses  at  a  fre¬ 
quency  of  10  megacycles  per  second  as  a  function  of  crystal  orientation 
and  of  temperature. 
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CRYSTAL  GROWING  AND  CRYSTALLINE  STRUCTURE  OE  BENZENE 


The  large  single  crystals,  necessary  in  order  to  follow  the  align¬ 
ment  and  cutting  procedures  developed,  were  grown  in  a  melt  of  reagent 
grade  benzene.  The  melt  was  contained  in  a  sealed  container  which  was 
maintained  at  a  temperature  of  approximately  3>°C.  After  degassing,  the 
bath  was  left  at  this  temperature  for  approximately  two  weeks  during 
which  time  the  crystals  formed. 

The  conglomerate  mass  of  crystals  was  then  broken  apart  and  large 
single  crystals  of  irregular  shape  were  obtained  and  stored  at  a 
temperature  of  about  -10°C  until  required. 

Crystalline  benzene  is  of  the  orthorhombic  class,  space  group 
Pbca,  the  unit  cell  having  the  dimensions  (2); 

a  =  7.292  A0 
b  =  9.U71  A° 

c.  =  6.71*2  A0  at  -195°C 

Optically  the  crystal  is  transparent  with  indices  of  refraction  (3) 

>  1.61* 

^  ^  =  o.ool* 

'6  -  ot  >  0.12 

The  crystals  were  thus  very  suitable  for  examination  by  optical 
methods. 

CRYSTAL  ALIGNMENT 

In  order  to  obtain  the  necessary  data  to  solve  the  strain  energy 
equation  and  thereby  compute  the  elastic  constants,  crystals  were 
required  which  were  accurately  aligned  along  each  of  the  principal 
axes  together  with  at  least  one  sample  which  could  be  aligned  along 
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some  arbitrary  axis. 

The  method  of  alignment  used  consisted  of  examining  the  raw  ir¬ 
regularly  shaped  crystal. s  under  a  conical  beam  of  polarized  light. 

For  this  purpose,  a  conoscope  was  constructed  in  the  form  shown 
schematically  (Fig.  1).  For  actual  examination  the  instrument  was 
located  in  a  refrigerated  compartment  maintained  at  a  temperature  of 
about  -10°C. 

In  the  raw  state,  since  there  seems  to  be  no  apparent  preferred 
direction  of  growth  for  benzene  crystals,  initial  examination  was  made 
with  the  crystal  immersed  in  a  glass  bottomed  tank  containing  a  mixture 
of  benzene  (CgHg)  and  methyl  salicylate  (CgHgOg) .  This  mixture  was  in 
the  proportions  of  Uj.0  parts  methyl  salicylate  to  300  parts  benzene  by 
volume.  The  refractive  index  of  this  mixture  is  approximately  1.32, 
close  to  that  of  benzene  at  -10°C.  The  actual  mixture  proportions  were 
found  to  be  not  critical.  Changes  occur  due  to  evaporation  and  the 
solubility  of  the  crystalline  benzene  in  the  mixture.  Fortunately 
the  mixture  is  close  to  the  solubility  limit  of  benzene  in  methyl 
salicylate,  so  that  loss  of  crystal  by  dissolution  is  small. 

The  tank  was  positioned  in  the  conoscope  so  that  the  bottom  of  the 
fluid  was  at  the  focal  point  of  the  convergent  lens,  with  the  divergent 
beam  spreading  up  through  the  fluid. 

In  this  beam  the  raw  crystal  was  examined  for  the  location  of  the 
optic  axes,  the  bisectrix  of  which  locates  the  b  axis  of  the  crystal 
(Fig.  8). 

Examination  was  carried  out  by  rotating  the  crystal,  held  in  a 
pair  of  tongs,  in  the  beam  to  determine  two  extinction  planes.  Since 
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the  normals  to  these  extinction  planes  form  an  orthogonal  system  with 
one  of  the  optic  axes,  this  axis  was  visible  at  90  degrees  to  both  these 
normal s . 

This  alignment  was  retained  by  transferring  the  crystal  to  a  special 
clamp  (Fig.  2),  designed  to  accommodate  a  range  of  shapes  and  sizes  of 
crystals  and  at  the  same  time  permit  a  more  accurate  alignment.  The 
clamp  consisted  of  a  short  length  (about  5cm)  of  10cm  diameter  aluminum 
tubing  into  which  were  mounted  in  gimbal  fashion  two  aluminum  rings.  On 
the  inside  of  the  inner  gimbal  were  secured  four  strips  of  spring  steel 
bent  into  the  center  forming  a  spring  collar  which  retains  the  crystal. 
Spring  pressure  was  just  sufficient  to  hold  the  crystal  firmly  in 
position. 

Mounted  in  this  clamp,  the  crystal  was  returned  to  the  bath  with 
the  base  of  the  clamp  forming  a  reference  face  resting  on  the  glass 
bottom  of  the  bath  and  was  thus  in  a  plane  normal  to  the  optic  axis  of 
the  conoscope. 

With  the  freedom  afforded  by  the  gimbals,  the  crystal  could  be 
rotated  with  respect  to  the  optic  axis  of  the  conoscope  to  bring  the 
bisectrix  into  coincidence  with  this  axis.  Since  the  crystal  was  held 
firmly  in  the  clamp,  the  interference  pattern  (Fig.  3)  could  be  clearly 
observed. 

The  alignment  was  preserved  by  removing  the  clamp  and  crystal  from 
the  bath  and,  using  the  base  of  the  clamp  as  reference,  melting  paral¬ 
lel  faces  on  the  crystal,  which  were  finally  covered  with  glass  slides. 

A  check  on  this  alignment  was  made  by  mounting  the  crystal  on  a 
Goniometer  attachment,  which  could  be  fitted  to  the  conoscope. 
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FIG.  3. 
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Small  corrections  to  the  alignment  could  be  made  by  melting  new 
faces  on  the  crystal,  the  corrections  being  computed  from  the  goni¬ 
ometer  angles  using  the  knowledge  of  the  refractive  indices  of  benzene. 

Alignment  to  the  other  principal  axes  was  accomplished  by  rotating 
the  goniometer  90  degrees  in  the  plane  of  the  optic  axes  (100),  or  in  a 
plane  normal  to  the  plane  of  the  optic  axes  (001) .  After  rotation,  new 
parallel  faces  are  melted  on  the  crystal  and  glass  slides  secured  to 
preserve  the  particular  orientation. 

Any  other  orientation  could  be  obtained  by  appropriate  rotation 
of  the  goniometer.  In  all  cases,  a  knowledge  of  the  direction  of  the 
optic  plane  was  preserved  by  appropriately  marking  the  glass  covering 
slides. 

MOUNTING  OF  CRYSTALS 

For  the  purpose  of  ultrasonic  pulsing,  it  was  necessary  to  mount 
the  crystals  in  such  a  manner  as  to  permit  electrical  connection  of  a 
radio  frequency  source  and  a  receiver  to  a  transducer,  which  was  in 
mechanical  contact  with  the  crystal.  At  the  same  time,  it  was 
necessary  to  ensure  that  the  crystal  orientation  was  maintained  and 
the  effects  of  sublimation  eliminated.  The  crystal  holder  (Fig.  U) 
was  therefore  designed  and  manufactured  to  meet  these  requirements. 

To  prepare  a  crystal  for  mounting  in  this  holder,  it  was  trans¬ 
ferred  to  a  glass  plate  and  frozen  in  place  using  the  preserved  face  as 
reference;  the  desired  axis  was  thus  maintained  normal  to  the  glass 
plate.  The  cylindrical  shape  required  to  fit  the  holder  was  then 
obtained  by  melting  the  periphery  of  the  crystal,  using  metal  cylinders 
maintained  at  room  temperature  to  melt  the  crystal .  When  the  desired 


8 


CRYSTAL 

BOLDER. 

SHOWING 

METHOD  OF 

MOUNTING 

AND  SEALING 

CRYSTAL . 

FI  G.  4. 

cylindrical  shape  was  obtained,  the  holder  was  slipped  over  the  sample 
and  frozen  in  place.  During  this  process,  the  sample  and  holder  were 
maintained  at  refrigerator  temperature. 

With  the  crystal  in  the  holder,  the  end  faces  were  finished 
parallel  to  the  holder  faces  by  polishing  on  fine  emery  cloth.  The 
end  faces  were  then  sealed  with  l/U  mil  clear  Mylar  for  the  free  end 
and  l/k  mil  aluminum  coated  Mylar  for  the  end  receiving  the  trans¬ 
ducer,  using  glycerine  as  a  sealing  agent.  The  end  seals  were  secured 
in  place  using  the  holder  end  plates. 

The  transducer  was  mounted  on  the  aluminum  coated  Mylar  using 
Nonaq  stopcock  grease.  These  bonding  materials  were  found  satisfactory 
for  the  transmission  of  both  shear  waves  and  longitudinal  waves  over 
the  temperature  range  for  which  data  was  obtained . 

The  transducers  employed  were  0.5  inch  diameter  quartz  crystal. 

X  cut  crystals  were  used  for  the  generation  of  longitudinal  waves 5 
Y  cut  for  shear  or  transverse  waves. 

PULSING  OF  CRYSTALS 

For  pulsing,  the  crystal  holder  was  mounted  in  the  pulsing 
assembly  (Fig.  5).  This  assembly  provided  a  secure  support  for  the 
holder  and  at  the  same  time  permitted  connection  of  the  transducer  to 
the  radio  frequency  source.  The  ground  connection  was  made  through 
the  body  of  the  assembly  and  the  crystal  holder,  the  other  connection 
being  made  through  the  spring  loaded  electrode.  The  method  of  mounting 
the  holder  allowed  adjustment  of  the  electrode  pressure  on  the  trans¬ 
ducer;  a  light  pressure  was  found  best  for  the  generation  of  longitudi¬ 
nal  waves  and  a  heavier  pressure  for  transverse  waves. 
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Radio  frequency  pulses  were  generated  by  an  Arenberg  pulsed 
oscillator,  Model  PG  65>0-C,  connected  by  low  impedance  coaxial  cable 
to  the  pulsing  assembly.  The  return  echoes  after  reflection  from  the 
free  end  of  the  crystal  were  amplified  by  an  Arenberg  wideband  ampli¬ 
fier  WA-600  and  displayed  on  a  Tektronix  Type  RMii^A  oscilloscope.  A 
Tektronix  Type  RM  181  time  mark  generator  was  used  for  control  of  pulse 
repetition  rate  and  for  external  sweep  control  of  the  oscilloscope. 
Figure  6  shows  a  schematic  arrangement  of  the  pulsing  system. 
MEASUREMENT  OF  VELOCITY 

To  obtain  accurate  measurements  of  the  velocity  of  the  ultrasonic 
waves  through  the  crystal,  a  modification  of  the  method  of  velocity 
measurement  described  in  (I4)  was  used.  Figure  7  shows  the  construction 
of  the  velocity  comparator  designed  and  built  for  this  modified  method . 
In  the  velocity  comparator,  pulses  generated  by  the  quartz  crystal 
transducer  immersed  in  the  liquid  medium  were  reflected  off  the  rigid 
boundary  and  received  at  the  transducer.  The  path  length  and  thus  the 
transit  time  of  the  pulse  was  made  variable  by  altering  the  position 
of  the  boundary  by  means  of  the  micrometer  screw. 

In  making  a  measurement,  the  pulse  trains  from  both  the  comparator 
and  the  crystal  were  superimposed  on  the  oscilloscope  and  the  path 
length  adjusted  to  bring  corresponding  echoes  into  coincidence.  Pre¬ 
cise  matching  could  be  obtained  by  observing  the  interference  effects 
as  the  echoes  were  brought  into  coincidence,  since  the  comparator  and 
benzene  transducers  were  pulsed  in  parallel. 

The  velocity  in  the  crystal  was  then  computed  from  the  known 
velocity  in  the  liquid  medium  by  the  relation: 
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where:  10  =  initial  micrometer  reading  corresponding  to  zero 
path  length. 

lc  *  length  of  path  through  crystal . 

1  =  micrometer  reading  at  coincidence. 

V]_  =  sound  velocity  in  the  liquid  medium. 

Vc  =  sound  velocity  in  the  crystal , 

The  liquid  medium  used  was  a  mixture  of  100  parts  distilled  water 
to  22  parts  absolute  ethyl  alcohol  by  volume.  This  mixture  transmits 
sound  waves  at  a  velocity  given  by  Willard  (5)  which  is  of  a  magnitude 
suitable  for  comparison  with  the  velocities  of  benzene  crystals  and  is 
substantially  independent  of  temperature  in  the  region  of  25°C. 
Temperature  control  of  the  device  was,  therefore,  unnecessary.  The 
temperature,  however,  was  monitored  and  corrections  made  for  the 
small  variations  which  occurred. 

Corrections  to  the  measured  length  of  the  crystal  were  necessary 
to  take  into  account  the  effect  of  thermal  contraction  of  benzene. 

These  corrections  were  made  from  the  data  on  the  thermal  expansion  of 
benzene  given  by  Cox  (2),  and  by  Euchen  and  Lindenberg  (6),  and  applied 
as  a  percentage  decrease  in  the  length  of  the  crystal. 

The  accuracy  of  the  velocity  measurements  was  dependent  on  two 
dominant  factors; 

(a)  Measurement  of  crystal  length;  this  length  was  measured  to 
.001  cm  and,  after  correction  for  thermal  contraction,  was  considered 
accurate  to  within  .01  cm  for  the  3  cm  path  length. 
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(b)  Pulse  matching  between  the  crystal  and  velocity  comparator; 
the  two  10  me  pulses  were  matched  to  within  one  half  of  one  wave  length, 
resulting  in  a  possible  error  of  .05yU  sec  in  12 yj  sec. 

The  possible  errors  introduced  by  other  factors  were  not  of 
significant  magnitude  in  comparison  to  these  two  factors.  Based  on  this 
estimate,  the  velocity  error  was  determined  to  be  less  than  1%. 

Measured  velocities,  using  the  same  crystal,  were  reproduced  to 
within  1%  over  the  temperature  range  250°K  to  170°K. 

TEMPERATURE  CONTROL  AND  MEASUREMENT 

Velocity  measurements  were  made  over  the  temperature  range  250°K 
to  170°K  by  cooling  the  crystal  mounted  in  the  pulsing  assembly  using 
liquid  air.  In  order  to  maintain  an  even  cooling  rate  the  sealed 
pulsing  assembly  was  inserted  in  a  silvered  Dewar  containing  liquid 
air.  Contact  with  the  liquid  air  was  made  through  a  heat  leak,  which 
permitted  control  of  the  cooling  rate. 

Temperature  was  measured  using  a  Copper-Con stantan  thermocouple 
placed  in  contact  with  the  crystal  holder. 

TEMPERATURE  RANGE 

The  upper  limit  250°K  was  predetermined  by  the  need  to  preserve 
the  crystal  from  sublimation  and  also  ensure  the  maintenance  of  a 
satisfactory  bond  between  the  transducer  and  the  crystal. 

Bond  failure  occurred  in  the  temperature  range  l50°K  to  170°K, 
thus  determining  the  lower  limit  of  170°K.  The  precise  nature  of  the 
failure  was  not  ascertained  but  it  was  believed  that  differential  con¬ 
traction  of  the  crystal  and  the  Mylar  end  seal  resulted  in  separation  at 
the  glycerine  bond  producing  a  void  between  the  transducer  and  the  crystal. 
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DETERMINATION  OF  ELASTIC  CONSTANTS 


For  a  crystal  of  orthorhombic  symmetry,  Group  3,  the  nine  inde¬ 
pendent  elastic  constants  may  be  arrayed  in  accordance  with  Cady  (l), 
as  shown: 
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Since  it  was  possible  to  align  and  mount  crystals  oriented  along 
each  of  the  principal  axes,  the  main  diagonal  terms  C-j^  of  the  matrix 
were  computed  directly  from  longitudinal  and  transverse  wave  velocity 
measurements  along  each  of  these  axes. 

The  crystal  orientation,  measured  velocity  and  associated  elastic 
constant  determined  by  this  method  were  as  given  in  Table  1. 


Table  1. 


DIRECTION 

OF 

WAVE  PROPAGATION 

DIRECTION 

OF 

WAVE  DISPLACEMENT 

VELOCITY 

MEASURED 

ELASTIC 

CONSTANT 

100 

vaa 

Gl'l  s  f>  Vaa2  | 

100 

010 

Vab 

C66  *  P  Vab2 

001 

Vac 

C$5  *  p  Vac2 

100 

Vba 

■  c66  "  p  Vba2 

010 

010 

Vbb  i 

;  C22  -  P  Vbb2 

1 

001 

} 

Vbc 

cuu  -  f  Vbc2 

ioo 

Vca 

c55  -  p  Vc  2 

001 

010 

Vcb 

QiU  =  P  Vcb2 

001 

Vcc 

C33  s  p  Vcc2 

Note:  For  axial  velocities,  the  first  letter  indicates  the  axis  of 

propagation,  second  letter  indicates  axis  of  wave  displacement. 
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To  determine  the  off-diagonal  terms  Cij  ,  the  transformation 


equations  (App.  I)  were  used.  Crystal  orientations  were  chosen  so  that 
each  off-diagonal  constant  could  be  computed  in  terms  of  a  measured 
velocity  and  certain  of  the  previously  determined  main  diagonal 
constants.  For  the  orthorhombic  crystal  the  desired  form  of  the 
equations  was  obtained  by  considering  a  rotation  in  an  axial  plane  to 
an  angle  of  U5  degrees  from  two  of  the  principal  axes  (Fig.  8).  This 
produced  three  different  orientations  giving  three  independent  sets  of 
equations,  each  of  which  could  be  solved  for  one  of  the  three  off- 
diagonal  terms. 

The  three  independent  sets  of  equations  were  obtained  by  sub¬ 
stituting  the  appropriate  values  for  the  direction  cosines  1,  m,  and 
n  in  the  equations  as  demonstrated  in  Appendix  I. 


Substitution  of  these  particular  values  reduced  two  of  the  Cj_j 
terras  to  zero.  The  remaining  terms  when  substituted  into  the  secular 
equation  for  irrotational  plane  waves  traveling  in  the  X/  direction 
gave  rise  to  the  simplified  determinant  solutions 


(Cii  -  />V2  )  (C-jj  -  pi2  )  -  Cij2  -  0. 


This  may  be  reduced  tos 


where  the  C  terms  are  defined  in  terms  of  the  elastic  constants  and 
the  direction  cosines  (App.  I).  This  results  in  an  equation  in  terms 
of  the  longitudinal  velocity  and  previously  determined  elastic  constants 
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The  crystal  orientation,  measured  velocity  and  associated  elastic 


constants  were  as  given  in  Table  2% 


Table  2 


Direction  of 

Wave  Propagation 

Velocity 

Measured 

ciJ 

Cii 

ckk 

110 

V110 

C12 

C11 

C22 

c66 

101 

V101 

C13 

C11 

c33 

Oil 

V011 

c23 

C22 

c33 

chh 

All  of  the  data  pertinent  to  the  computations  and  the  actual 
values  of  the  elastic  constants  have  been  presented  in  both  graphical 
and  tabular  form,  page  21-26. 

ERROR  IN  ELASTIC  CONSTANTS 

The  sources  of  error  in  the  determination  of  the  elastic  constants 
from  the  measured  velocities  could  be  divided  into  two  categories  -  the 
error  in  crystal  alignment,  and  the  error  in  computing  the  elastic 
constants  from  the  velocity  data. 

The  error  in  crystal  alignment  depended  not  only  on  the  accuracy  of 
the  optical  alignment,  but  also  on  the  method  of  preserving  this  align¬ 
ment  during  the  mounting  procedure.  Since  for  all  alignments  other  than 
that  along  the  b  axis,  the  alignment  of  the  crystal  could  not  be 
rechecked  after  mounting,  no  correction  could  be  applied. 

Another  source  of  alignment  error  was  introduced  in  the  alignment 
of  the  shear  transducer  to  produce  shear  waves  along  the  axes.  To 
eliminate  this  error,  only  longitudinal  wave  velocities  were  used  to 
determine  the  off-diagonal  elastic  constants. 


20 


TABLE  OF  MEASURED  VELOCITIES  (meters  per  second) 


T°K 

Vaa 

Vab 

v 

ac 

vba 

vbb 

vbc 

vca 

vcb 

Vcc 

V110 

V101 

v011 

250 

231*1* 

11*12 

1888 

1362 

21*87 

1200 

1832 

1193 

2l*0l* 

2572 

2879 

2530 

21*0 

2381 

11*69 

19U0 

11*19 

2567 

1216 

1897 

1207 

21*55 

261*1* 

291*6 

2611 

230 

21*26 

1519 

1981* 

11*79 

2628 

121*1 

19U7 

1227 

2505 

2727 

3009 

2675 

220 

21*87 

1569 

2038 

1532 

2708 

1258 

1990 

121*7 

2550 

2796 

3069 

2737 

210 

2536 

1617 

2091 

1576 

2767 

1275 

2031 

1261* 

2593 

2865 

3130 

2793 

200 

2580 

1659 

2139 

1618 

2821 

1295 

2072 

1281 

2632 

2921 

3185 

281*9 

190 

2629 

1702 

2188 

1660 

2862 

1311* 

2113 

1300 

2675 

2975 

321*1 

2899 

180 

2677 

1735 

2237 

1693 

2905 

1336 

211*7 

1317 

2713 

3026 

3291 

291*3 

170 

2725 

1773 

2281* 

1731 

2951* 

1355 

2181 

1333 

271*7 

3071 

3350 

2988 

VELOCITY  (METERS  PER  SECOND 


AXIAL  VELOCITIES 


velocity  (Meters  per  second  x  io_ 


LONGITUDINAL  VELOCITIES 


ELASTIC  CONSTANTS  OF  BENZENE  (C6H6) 

x  10"^  dynes/cm.^ 


Temperature 


°K 

C11 

C22 

c33 

caa 

c55 

c66 

c12 

C23 

c13 

250 

5.83 

6.56 

6.1U 

1.53 

3.78 

1.97 

3.90 

3.52 

a. 01 

2J4.O 

6.01 

6.99 

6.39 

1.57 

3. 99 

2.ia 

a. 06 

3.57 

a. 22 

230 

6.2U 

7.33 

6.66 

1.63 

a. 17 

2.32 

a. 3a 

3.67 

a.ai 

220 

6.56 

7.78 

6.90 

1.68 

a.ai 

2.a9 

a.aa 

3.68 

a.a3 

210 

6.82 

8.12 

7.13 

1.72 

a. 6a 

2.63 

a. 67 

3.70 

a. 53 

200 

7.06 

8.UU 

7.35 

1.78 

a. 85 

2.78 

a. 76 

3-76 

a. 61 

190 

7.33 

8.69 

7.59 

1.83 

5.08 

2,92 

a. 92 

3.89 

a.  66 

180 

7.60 

8.95 

7.81 

1.89 

5.31 

3.oa 

5.07 

3.9a 

a.  66 

170 

7.88 

9.26 

8.01 

1.95 

5.53 

3.18 

5.08 

3.85 

a. 80 

2k 


ELASTIC  CONSTANTS  (Cjj) 


i 


ELASTIC  CONSTANTS  (  C :  ) 

*  j 


For  the  crystals  not  aligned  on  a  major  axis,  a  further  source  of 
error  was  introduced  by  the  anisotropic  thermal  contraction,  which 
produced  a  continual  minor  change  of  alignment  throughout  the  tempera¬ 
ture  range.  Because  of  other  larger  sources  of  error  present  in  the 
determination  of  the  off-diagonal  constants,  no  correction  was  made 
for  this  error. 

The  specific  volume  of  benzene  was  considered  to  be  0.9k2%  cc/gm 
for  the  temperature  range  25>0°K  to  170°K,  based  on  data  summarized  in 
Andrew  and  Fades  (8).  Possible  errors  in  this  data  were  considered 
negligible  in  comparison  with  the  other  errors  present. 

Since  it  was  not  possible  to  measure  any  alignment  errors  present, 
an  estimate  of  the  accuracy  was  obtained  by  comparison  of  velocities 
which  determine  the  same  elastic  constant.  For  the  constants  determined 
by  longitudinal  axial  velocities,  the  error  was  estimated  at  less  than 
2%,  based  on  velocities  obtained  from  different  crystal  samples.  For 
the  constants  determined  by  shear  velocities,  the  comparative 
velocities  obtained  also  differed  by  no  greater  than  2$. 

The  values  of  the  diagonal  elastic  constants  (C^)  determined  from 
axial  velocities  were  required  in  the  solution  of  the  off-diagonal 
elastic  constants  (Cij)$  thus,  any  error  present  in  these  values  was 
unavoidably  included  in  the  off-diagonal  values.  This  resulted  in  a 
possible  overall  error  in  the  off-diagonal  elastic  constants  of  up  to 


10% . 

Table  of  Estimated  Error  in  Elastic  Constants 


Elastic  Constants 

Velocities  Used 

Error 

cll>  C22,  C33 
cIiH’  c5!?>  c66 

Cl2»  C13,  C23 

Longitudinal,  Axial 

Shear,  Axial 

Longitudinal,  Non-axial 

less  than  2% 
approximately  2% 
up  to  10$ 
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MEASUREMENT  OF  ABSORPTION 


Although  measurement  of  absorption  was  not  the  primary  objective, 
the  method  used  to  pulse  crystals  to  obtain  velocity  data  was  readily 

adaptable  to  providing  absorption  data. 

Absorption  measurements  were,  therefore,  made  in  conjunction  with 
velocity  measurements  by  comparing  the  amplitude  of  successive  pulse 
echoes.  The  comparison  was  made  by  matching  the  leading  echo  to  the 
output  of  the  time  marker  generator  displayed  on  the  second  channel  of 
the  oscilloscope.  The  time  marker  generator  output  was  fed  through  an 
attenuator  which  permitted  an  adjustment  of  amplitude  in  steps  of  one 
decibel . 

In  measuring  the  attenuation  of  successive  echoes,  the  time  marker 
was  matched  to  the  leading  echo  and  the  attenuator  then  adjusted  until 
the  time  marker  amplitude  was  the  same  as  that  of  the  succeeding  echo. 
The  attenuation  between  echoes  was,  therefore,  given  directly  as  the 
attenuator  setting  in  decibels.  These  measurements  were  repeated  for 
as  many  echoes  as  appeared  in  a  given  pulse  train. 

All  absorption  measurements  were  made  at  a  frequency  of  10  mega¬ 
cycles  per  second  using  0,50  inch  diameter  quartz  crystal  transducers. 
The  accuracy  of  the  measurement  was  extended  to  within  one  half  of  one 
decibel  by  logarithmic  interpolation  of  the  superimposed  waveforms. 
DETERMINATION  OF  ABSORPTION 

The  average  of  the  attenuator  settings  was  then  converted  to 
absorption  in  nepers  per  centimeter  by  a  calculation  based  on  the 
measured  length  of  the  crystal: 

Absorption  ( neper  s/cm)  =  /f.32  x  8.68  nePers 

2  x  lc  (cm)  db 
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The  ultrasonic  absorption  computed  from  the  measurements  of 
pulse  attenuation  has  been  presented  graphically  on  page  30  and  in 
tabular  form  on  page  31. 

All  of  this  data  was  based  on  the  attenuation  between  first  and 
second  echoes  only,  thus  involving  a  total  path  length  of  six  centi¬ 
meters  in  the  crystal.  For  this  length  of  path,  the  effects  of 
diffraction  were  small  and  within  the  errors  of  the  actual  measure¬ 
ments,  thus  no  correction  was  made  for  attenuation  due  to  diffraction. 

The  absorption  computed  varied  from  greater  than  0.6  nepers  per 
cm.  to  less  than  0.1  nepers  per  cm.  the  variation  being  a  definite 
function  of  crystal  alignment  in  relation  to  the  direction  of  ultra¬ 
sonic  wave  propagation.  The  absorption  for  all  waves  propagating 
along  the  c  axis  was  significantly  less  than  for  waves  propagating 
perpendicular  to  this  direction. 

Some  variation  of  absorption  with  temperature  was  observed  for 
certain  of  the  crystal  directions.  In  general  this  variation  showed 
a  decrease  in  absorption  with  a  decrease  in  temperature . 
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ULTRASONIC  ABSORPTION 


ULTRASONIC  ABSORPTION  IN  BENZENE 


Direction 

of 

Propagation 

100 

010 

110 

101 

010 

100 

001 

Oil 

001 


Temperature 

Range 

°K 

250  -  170 

250  -  230 
230  -  170 

250  -  230 
230  -  170 

2U0  -  225 
225  -  210 
210  -  195 
195  -  170 

250  -  230 
230  -  210 
210  -  200 
200  -  170 

250  -  21*0 
21*0  -  220 
220  -  170 

250  -  21*0 
21*0  -  170 

250  -  215 
215  -  205 
205  -  170 

250  -  210 
210  -  190 
190  -  170 


Attenuation 

db/cm. 

5.3 

3.3 
U.3 

6.0 

3.5 

a.7 
a. 2 

3.5 

3.0 

a. 6 
a. 2 

3.7 
3.0 

3.3 
3.0 

2.7 

2.0 

1.7 

1.7 

1.5 

1.3 

2.3 

1.3 
0.7 


Absorption 
nepers 
per  cm. 

0.61 

0.38 

0.a9 

0.70 

o.ao 

0.5a 

o.as 

o.ao 

0.35 

0.53 

o.as 

0.1*3 

0.35 

0.38 

0.35 

0.31 

0.23 

0.19 

0.20 

0.17 

0.15 

0.26 

0.15 

0.08 


Mode 

Shear  (001) 
Shear  (001) 

Long. 

Long. 

Long. 

Long. 

Long. 

Long. 

Shear  (010) 
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APPKNDIX  I 


The  method  developed  by  J,  R.  Neighbours  (7)  for  the  determination 
of  elastic  constants  is  based  on  velocity  measurements  of  irrotational 
plane  waves  propagated  along  a  specified  reference  axis  (X/),  where  the 
reference  axes  are  oriented  so  that  the  Y'  axis  lies  in  the  001  (XY) 
plane  of  the  crystal  . 

Using  the  measured  velocity  data,  the  elastic  constants  are 
determined  from  the  solution  of  the  secular  equation: 

(Cji  -/’v2)  °16  c15 

cl6  (c66  -f’'’2)  C#  -  0 

C15  <=56  (c 

where  the  primed  constants  are  defined  in  terms  of  the  crystal  elastic 
constants  by: 

Cqi  -  Z  4  2£.^m^C^2  4  2d2n'LC^  4  m^C22  +  2m2n‘~C23  4  n^C^^ 

4  Iim2n2C^  4  U^2n2C^  4  U£2m2C66 

=  (2m2a-2C^q  -  2X2m2a“2C]_2  i  C  2m2a“2C22  +  X  2n2a"2C^ 

4  m2n2a"2Cc£  4  (/2  -m2)2a~2C^ 

C55  =  X^n'  a~2C]_2  “  2j  2m2n2a~2Cq2  ~  2L2n2C-^^  4  mXn2a"2C22  -  2m2n2C23 

4  n2a2C33  4  m2a"2(a2-n2)2C[t^  4  X2a~2(a2-n2)2C^  4  lj^/2m2n2a“2C66 
0x5  =  ■'X^na-^Cxx  “  2  X2m2na”^Cx2  “  C2na“-*-(a2-n2)Cx3  -  nXna'-k^ 

4  ra2na~l(a2-n2)C23  4  0^033  4  2m2na“-*-(a2-n2)C[lj1  4  2  X  2na”^(a2-n2)C^^ 
-  h  ( 2m2na“-*-C66 

cl6  =  -X^a_1C3x  +  /ma“1(-62-.m2)C12  -  Xmn^a^C^  4  X  nXa_1C22 

4  Cmn2a"-*-C23  ^  2  f  mn2a~‘-*-C)|).  -  X  mn2a“^-C^  4  2  1  ma"-*-(X  2-m2)C^^ 

C56  =  ,t 3mna“2Cxx  -  Xmna”2(.X  2-m2)C-L2  -  ^mnC-^3  -  ^m^na”2C22 

4  (mnC23  +  •^mna“2(a2-n2)Cj1j1  -  mna-2(a2«n2)Cc£  -2/mna“2(i  2-m2)C6£ 

inhere  a2  =  C 2  4  m2 
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and  the  C,  m,  and  n  terns  are  the  direction  cosines  relating  the  X 
direction  to  the  crystallographic  axes  (a,b,c). 

For  a  wave  propagating  along  a  crystallographic  axis,  the  direction 
cosines  become  one  or  zero,  which  results  ins 

Ci  j  -  o 
cii  =  cii 

so  that  the  solution  is  in  terms  of 

pi2  s  Cii  o 

Example :  f or  4  -  1  Cu  =  Cn  *  ^V2aa 

m  =  n  -  0  c£  =  Ctf  -  pl2ab 

c66  =  C66  *  ^V2ac 


For  a  wave  propagating  in  an  axial  plane  at  hS  degrees  to  two 
crystallographic  axes,  the  direction  cosines  become  zero  or  V2/2* 
For  the  110  case,  where  the  wave  is  propagating  at  l|f?  degrees  to  the 
a  and  b  axes,  and  perpendicular  to  the  c  axiss 


-C  =  m  -  42/2 

cn  2  Kci:l  t  C22  4  ^c66  +  2^12) 

P 

IB 

O 

c66  "  4 (C11  4  c22  -2C12) 

a  s  1 

c55  "  KcItU  4  c5^) 

C16  “  KC22  “cll) 

Cl5  s  C56  25  0 


thus  the  determinant  becomes: 

C11  ~  cl6  0 

C16  (c66  "  P^2)  0 

0  0  (C#  -  pi2) 

or  (Cjj  -^V2)  [(Ch  -^V2)(C6g  -^>V2)  -  C^2]  =  0 


=  0 
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which  has  solutions: 


(a)  pV2  =  0^  »  i(CW,  t  C35) 

(b)  Pv2  »  Men  1  c66)  ±  [Ken  -  c6'6)2  1  c^2] 

The  solution  for  the  longitudinal  velocity  (V^  )  is: 


p v  ,  2  =  l(CU  -  c66)  ^  [4(^11  -  o{6)2  i  <tf]  2 

=  4^11  +  ^22  ^  ^c66 ^  i  [(C12  +  c66^2  +  ^(’-'22  “ 


Separating  the  off-diagonal  elastic  constant 


'12 


(cll  t  C22  +  2C66)]  -  i(C22  -  Cn)2/ 


c66 


results  in  a  solution  for  the  off-diagonal  constant  in  terms  of  the 
longitudinal  velocity  and  the  previously  determined  diagonal  constants. 
The  solutions  for  the  other  off-diagonal  constants  are  analagous: 


for  the  101  Case  C 


13 


for  the  Oil  Case  C 
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2  ~  ?(cll  +  c33  +  2C55)]~  4(c33~cn)  j"  2  -  C33 
2ypV£2  -  |(C22  +  t  ^33)]"  4(<333”C22)2^  2  _ 


It  may  be  noted  that  since  the  solution  is  of  the  form: 


(C^4j  -f-  Ckk)  -  a  known  positive  constant  (k2), 
this  allows  for  the  possibility  of  two  solutions: 

cij  ~  K  -  Ckk 
or  Ci2  -  -(k  +  Ckk) 

The  choice  between  solutions  was  made  on  the  basis  of  the  stability 
requirements  for  a  crystal  of  orthorhombic  symmetry.  From  Neighbours  (9) > 
this  requires  that  the  principal  minors  of  the  matrix  must  all  be 
positive.  The  matrix  was  therefore  expanded  into  terms  of  the  principal 


minors,  and  the  following  four  inequalities  obtained? 


cllc22  >  c122 
cllc33  >  C132 
c22c33  >  c232 

C11^C22C33"C2 32^  +  2C12C23C13  ^  c122g33  f  C132c22 
These  inequalities  were  then  used  to  show  that  the  solution  giving  a 

large  negative  value  was  not  acceptable. 
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